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r~i , Abstract 

. 

An algebraic method of constructing the confluent Natanzon potentials en- 
^ ' dowed with position-dependent mass is presented. This is possible by identify- 

ing the scaling resolvent operator (Green's function) to nonrelativistic position- 
dependent mass Schrodinger equation. 
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£4 1 Introduction 
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The (confluent) Natanzon potentials [1,2] are known nowadays as a general class of 
solvable potentials related to the (confluent) hypergeometric functions and cover all 
well-known potentials in quantum mechanics for which an exact solutions can be 



>< 



found [3]. The (Radial) harmonic oscillator potential [4], the Coulomb potential [4] 



and the Morse potential [5] may all be considered as particular cases of the confluent 
hypergeometric potentials, while the hypergeometric potentials consist of the Poschl- 
Teller potential [6], Eckart potential [7], Rosen-Morse potential [8], etc. 

However, different series of papers have been devoted in recent years to the study of 
the dynamical symmetries [9-17]. It was proved that the bound and scattering states of 
these potentials are related to the representation of the so(2, 1) [~ su (1, 1) ~ si (2, R)] 
Lie algebras; these algebras play the role of dynamical algebras. They are used as a 
spectrum generating algebra [9-13] and as a potential algebra [14-17] in order to obtain 
the analytic solutions. 
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We have succeeded in generating the Natanzon hypergeometric potentials by mak- 
ing use of conformal mappings within the framework of position-dependent mass 
(PDM) [18]. In this paper, we will explain the form in which the so(2, 1) Lie al- 
gebra can be enlarged as wide as possible in order to generate the confluent Natanzon 
hypergeometric potentials endowed with PDM. This is accomplished by scaling the 
resolvent operator (Green's function) [19], in such a way that can be related to non- 
relativistic PDM Schrodinger equation [20-24]. This explicit connection between the 
resolvent operator transformations to PDM quantum systems. 

The organization of the present paper is as follows. We introduce in section 2 a 
specific differential realization of the so (2, 1) Lie algebra. Then applying a scaling in 
the resolvent operator, through a scale function that preserves the wavefunctions, will 
map it into the wave equation endowed with position-dependent mass. This is possible 
if the scale function is identified to the mass distribution. This correspondence gives, 
in section 3, the (effective) potential expression for confluent Natanzon potentials in 
the framework of position-dependent mass. The last section is kept for our concluding 
remarks and an appendix was added in order to point out the correct use of the scalar 
product under this scaling. 



2 Resolvent operator transformations and 50 (2, 1) 
algebra 

We are going to build a realization of the so (2, 1) Lie algebra according to [12]. Let 
us start from a couple of operators a and b satisfying the commutation relation 



[a,b] = l, 

and defining the following realization for the generators J k (k = 0, 1, 2) 



Jo = 4 



Jl = 4 



fl 2+ 4c+t b- 2 + b 



a 2 + 4c + - 6 



J 2 = __( 2 &a + l), 



(1) 

(2.1) 
(2.2) 
(2.3) 



where c are eigenvalues of the Casimir's operator of the algebra. The relations (2) 
satisfy the commutation relations 



[Jo,Ji]=U2, [J2,Jo]=iJi, [Ji, J2] = -iJo- 



(3) 



It might be of interest to consider the operators a and b in terms of a variable x as 
a = -£: + K (x) and b = x. Performing the point canonical transformation x = a/£ (u), 
the operators a and b become 



a = 



u 



?(u) du 2 I 



(4) 
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where here and below the prime stands for derivatives in variable u. In terms of this 
realization, the generators (2) become [12] 

_ e{u)du* 2 e( U ) + 4 ^( u ) £w 4 ' 

1 " rw^ 2 2 e /3 n 4 r(«) ew 4 ' 1 ' 

The Green's function for a wide class of potentials is given following expression 

(H-E)G(u,u)=6(u-u), (6) 

and due to the form of the generators in (5), it is easy to give the resolvent operator 
[19] 



= qo + ^2?kJk, (7) 



fc=0 

where go an d Pk (k — 0, 1, 2) are arbitrary parameters that will be fixed below. 

The realization in (7) may be implemented by scaling the resolvent operator on 
the generators (2) in such a way to explicit connection between the resolvent operator 
(7), and then the spectrum generating algebra method, to PDM Schrodinger equation. 
Under such transformations, the operator Q (E) becomes 

Qo + ^PkJk P(u,u) 




= P( u ,u)^-(H-E (8) 

where P (u, u) = P± (u) P2 (u) is called two-point scale function, and Tk are generators 
which are determined through identification. 

Under this change and by identifying both sides of (8) in terms of parameters q 
and pk, the generators T k read as 

° £ 2 d U 2 fPdu " 2 » o ,/3 + , ,/4 + , + 



( P" 




3 « 
4 r 


c 

+ r 


e 




2 £' 3 


4 




1 re 


3 re 

+ 4 r 


c 


e 




2f 


4 



1 fdv? fPdu fP 2r 4f e 4' 1 j 

= -^"Vp'ar (9 - 3) 
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where P here is referred to the function Pi (u) due to differentiation with respect to 
the variable u. 

By introducing the eigenfunctions 



<P{u) 



(10) 



the operator Vt (E) becomes after some straightforward algebras 



n(E) = p(u,u)^ 



' 2p+ ^& 2 " I { 4p+ ( 2 + ^ ~ 2 ^ 1 + 2mi 



p, t pn t2 

r (8-^ + 4-20 + 2— ^ „ p ^ 



I /-2/-/I t-fil oi f2<-/l2 



d_ 

du 

2 /-I" 



*/2 



2 e /4 
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+2qot + yf - 2* P2 (J^ + 2£ - ~- ^ 
where p± = po ± p\ . 

Without any loss of generality, let us take the appropriate-parameter choices on 
p± and p 2 such as 



3£"£ 2 



(11) 



- , a - & -n 

P+=Po+Pi=2' P-=Po-Pi=2> ^2 = 0, 



(12) 



then the operator f2 (i?) in (11) becomes 



ft(£) = P(u,u)-^[H-E 



= P(u,u)-^ 



-a 



£' 2 du 2 



"77 4 + 2 



p_i_ A ?z\ d 



du 

5/ tf 2 e // 



+2%^ + ^ 2 - 4a ^^ - 2a + ac - ^ + 4a p-^t + 6a ^- 
21a i 2 i" 2 3aZ 2 g"~ 



4 e 



2 e 



/3 



(13) 



On the other hand, in the case of varying mass which will be denoted by M (it) = 
m m(u), the Hamiltonian introduced by von Roos [20] reads as 



H v r = ^ [mPprripmP + m p pm e pm v ) + V c s (u) , 



(14) 



where r], e, and p are three parameters which obey to the restriction rj + e + p = —1 in 
order to grant the classical limit. Here p (= —i-^) is a momentum with /i = mo = 1, 
and m (u) is dimensionless-real valued mass. 
Applying the eigenfunctions [18,23] 



^ (w) = m (w) ^ (it) 



2,,/U)4M0(u) 



r («) 



(15) 
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the resolvent operator of H v r reads [18] 

H vR -h = -775— - 77 I 4 + —— —)— + TJ2\ 6 ^ +— TFT 



£' 2 du 2 £' V mg i' 2 J du £ 2 \ £' £ 

m £ \ ££ 2 m' m 

£ 2 

+2m|__(V eff (x)-^). (i 6 ) 
Then comparing (13) to (16) in their differential terms, we obtain 



P (u,u) = \J m (u) m (uj, (17) 

leading to defining the one-point scale function P (u, u) = P (u) — m (it) as it = u. 
By substituting (17) into the remainder term leads, after long calculations, to the 
following general expression for the potential 

v K <«>1 " E = 2^) fw ( 2 * ? M + ^ M + c ) ~ 5STO0 {? W '" >s ' (18) 

and which corresponds to the specific-parameter choice a = 1. Here, {£ (it) , -u} s is 
the Schwarz derivative of £ (it) with respect to it defined by 

{*(«) «> =^P-^1 (19) 

and 1/ [£ («)] = V eS [£ (u)] - V&' e) (u), with 

Mt) _ (l + 2r ? ) 2 + 4e(l + ? ? ) m' 2 (u) _ e m" (u) 

Vm W ~ 8 w?{u) 4m 2 ( M )' { ' 

which is attributed to the dependence of the mass m on u. 



3 Generating confluent Natanzon potentials 

Since we have to get a constant term in left-hand side of (18), we must divide the 
right-hand side into two parts: the former represents a constant term and which is 
identified to the energy E, while the latter corresponds to the potential V [£ (it)]. This 
condition amounts to setting up differential equation for £ (it) and satisfying 

i' 2 (u) A 2 £ 2 (u) + Ai£ (u) + Aq 

2m (u) 4£ 2 ( W ) ' 1 J 

where Aj (i = 0, 1, 2) are arbitrary parameters. Inserting (21) into (18), we get 

v K («)] - * = ^>)±MM±ig _ 1 {e (tt) , M}s . (22) 

1 1 ;J A 2 £ 2 (it) + A : £ (u) + A 4m (it) s ^ 
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However in order to get a constant term in right-hand side of (22), it is useful to 
make some particularizations on the parameters f3, qo, and c. As a consequence, it 
will turn out to be useful to express the resolvent operator — E^j in (13) as the 
exponential function of the generator T 2 , since T 2 is essentially a generator of scale 
transformations [25]. 

To this end, taking into account (12), the resolvent operator (8) reads 

(h-E)iP (u) = [2 (f3 + 1) T - 2 (/? - 1) 7\ + 8g ] ^ («) = 0. (23) 

It follows from the identity [26] 

oo 1 



J0T 2 



(n) 
fc=0 

T cosh - 7\ sinh 0, 



(24) 



where 



[i0T 2 ,T o ] (n) = [ieT 2 ,[ieT 2 ,T ] in _ 1} 

= i9T 2 , i9T 2 ,[i9T 2 ,T ] {n _ 2) 



that the right-hand side of (23) becomes 

%ifj (u) = e im T e~ i5T2 $ (u) 
where 9 and 5 are expressed in terms of (3 as 



2?o 

VP 



9 = j-^-, and 5 = ^\n(3, 



(25) 



(26) 



(27) 



and are allowed to depend on n and c. 

Since J is the only compact generator with eigenvalues ( J ) = jo> then the allowed 
values of j are related to n and c through [12] 

Jo = n+i(l + v / TT4^), (28) 
where n = 0, 1, 2, . . . Combining (28) to (26), we get the basic equation 

VP 



0. 



(29) 



It is then obvious that the simplest set of parameters associated with f3, q and c 
and compatible with (29) is 



[3 = a 2 , q = - 



a ( b 



2 \2 



+ n 



b (b 



2 \2 



- 1 



(30) 
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where the equation (29) can be obtained using the equalities (30) by the straight 
substitution. 

Equation (22) becomes 

V K («)] -E = M - 2« (2n + » Hu) ± b £ - 2) _ 1 

which in particular case requires that coefficients in numerator are linear with those 
of denominator with respect to E following 

a 2 = -\ 2 E + ap, (32.1) 
2a(2n + 6) = X ± E - a qo , (32.2) 
6(6-2) = -\ E + a c , (32.3) 

leading to express the potential V (u) into 

^ (m)] = — i^y {e(M) ' M}s ' (33) 

with R (f ) = A 2 £ 2 + Aif + A . For convenience, we introduce the auxiliary function 
from (21) as [12,18] 

6[ ^ M)] = 2^)=^M' (34) 
therefore the Schwarz derivative-term can be expressed in terms of £ (it) 

(«) , u} s = -- — + ^^~U m («) , (35) 

where the dot refers to the derivative with respect to £ and U m (u) = Jj^^Vrr — | m 2 /^ - 
Inserting (35) into (33), taking into account the first and the second derivative, we 

get 

T/ / x _ ^ + ggo£ + ^ + 1 , 4A 2^ 2 + Aig 5 c2 (4A 2 e + AQ 2 

S(g + ^ 2 (0 i e i? 3 (0 ' (36J 

where a effective-mass potential is now given 

W#' e) (u) = Vt e) (u) + U m (u) 

4(l + 27]) 2 + 16e(l+?]) + 5m /2 (n) 2e + lm"(«) 

32 m 3 (u) 8 m 2 (u) ' ( ^ 

Now consider simultaneous addition and subtraction of A 2 — A A 2 , this procedure 
allows us to get R(£) in the numerator of the third term. We obtain after some 
straightforward algebras 

1 {u) = 5(0 + ~iP(£) I* *Hi) • (38) 

where A = A 2 — 4AoA 2 . We recognize in (38) the general expression of the confluent 
Natanzon potentials. Note that the potential deduced is completely independent of 
ambiguity parameters rj, e, and p. However, the potential differs substantially from 
that usually known in the literature [11,12]. 
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3.1 Energy eigenvalues E n 

The corresponding energy eigenvalues of the confluent Natanzon potentials can be de- 
termined from (32). In fact, proceeding to squaring (32. a) and (32. c) and by inserting 
it into (32. b), we see then that 

-cr Q o + Ai-E 1 , 



2n + 1 = ^ T^V^ -Vl + vc- AqK, (39) 

2 V a l3 ~ *2&n 

is an equality of the fourth degree in E n and agree with that obtained by [11,12]. It is 
obvious that parameters a and b depend on quantum number n. 

3.2 Wavefunctions ip n [£ {u)\ 

The two sets of operators Tj and % are obviously different but have the same effects 
on the wavefunctions. Using (24) and (26), it turns out that the parameter = a n 
plays the role of scaling parameter in such a way that £ (u) can be mapped into a„£ (u), 
but does not affect the general form of wavefunctions. 

The wavefunctions belong to the Hilbert space, not with respect to the usual scalar 
product of quantum mechanics but with respect to the new scalar product (see the 
appendix) 



/+oo 
fa ( u ) i>] ( u "> du 
— oo 



/ 



oo 
+oo 

— oo 



Xi ( u ) fp^yXj («) du , (4°) 



where ip^u) = y=pXi ( u ) (« = 1,2,...) and W (u) = = ^ represents the 

weight factor. 

Then, and following [9] , the assumption of an unnormalized wavefunction ip n [£ {u)\ 
gets transformed as 



^„K(«)] ~ ^^e n '\u)e-^/\F 1 {-n-b n -a n au)) 

= m 1 / 4 (u) R 1 ^ (u) C (K - l)/2 (u) e-««)/ 2 ^ (-n; 6 n ; a„£ («)) , (41) 

where the identity (34) is used. The wavefunction (41) has an extra term a/ m (u) 
compared with that for the constant mass (see the appendix). This achieves the 
bound-state problem for the confluent Natanzon potentials. 



4 Conclusion 

We have investigated another method to generate the confluent Natanzon potentials 
and its bound-states. The method consists to enlarge the so (2, 1) Lie group genera- 
tors in order to relate a resolvent operator transformations to the position-dependent 
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mass Schrodinger equation. We have shown that if a particular correspondence be- 
tween the scale function and mass distribution exists, i.e. P (u, u) = ^Jm(u)m (u), 
then there is a connection procedure between the resolvent operator transformations 
and position-dependent mass Schrodinger equation leading to generate the confluent 
Natanzon potentials and their bound-states. We have also pointed out the correct use 
of the W (w)-scalar product when dealing with position-dependent mass. 



5 Appendix: The correct scalar product 



As outlined in the introduction, this appendix was added in order to bring some 
mathematical details were not approached in subsection 3.2 relating to the correct use 
of the scalar product under the P (-u)-scaling. 

The time-dependent Schrodinger equation for the functions ip (u, t) and ip* (u, t) 
and associated to the Hamiltonian (14) read, respectively, as [27] 



.dip (u, t) 
.dip*(u,t) 

1 Ft 



d 



d 



du 2m (u) du 
d 1 d 
du 2m (u) du 



i>(u,t) + V ef{ (u)4){u,t), (A.l) 
V>* (u, t) + V eS (u) i>* (u, t) , (A.2) 



where V e s (u) G K. 

Let tp 1 (u,t) and il)* 2 {u,i) be the solutions of (A.l) and (A.2) respectively. By 
multiplying (A.l) by W (u) ip*, (u, t) and (A.2) by W (u) ip 1 (u, t) on the left-hand side, 
we end up on subtracting both equations to the equality 



ij^i (u,t)W( U ) i>*2 («,*)] : 

d 1 



W(u) 



ip* 2 (u,t) 



du 2m (u) du 



d did 

i>\ ( u ,t) ~4>i (u,t) 



du 2m (u) du 



1p2 (U,t) 



, (A.3) 



due to the fact that W (u) is independent of time t. Then (A.3) can be recast into the 
differential form as 



(u,t)w( U )r 2 («,*)] 

= -»W(u)l<W(u) 



du 



du 



, (A.4) 



with a specific choice of the weight function W (u) chosen as 
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W(u) = 



NP(u) Nm(u)' 



for any N E N - {0} . 



(A.5) 



Let us introduce the wavefunctions Xi ( u i t) (i — 1, 2, . . .) related to the eigenfunc- 
tions (u, i) by 

Xi (u,t) = y/W{u)^ (u,t), (A.6) 
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then (A. 4) can be rewritten as 
■ d ( ( v 1 A N d ( dxi (u,t) d X * 2 {u,t) 

Xl K t) 777^X2 ( U ^) ) = -1777- \X2 (U, t) Xl K t) 



(A.7) 

By integrating (A.7) with respect to the variable u, we get 



d f +co 1 

»L «' M fi5« !M '" = * (A ' 8) 



where, as is excepted for bound-state wavef unctions, Xi ( u i t) ~ ¥ an d X2 ( M ) t) —> 
as ti — > ±00. It becomes clear that these wavefunctions are not orthonormal as they 
stand; rather they are orthonormal with respect to the weight factor w \-. defined in 
(A. 5). Dealing only with spatial wavefunctions, i.e. Xi( u )i (A. 8) leads to the new 
scalar product 

/+00 ^ 
Xi (u) ^T(u) X * 2 ^ du ' ^ A ' 9 ^ 

which is more general than that given usually in quantum mechanics. Through (A.9), 
taken into account (A. 5) and N — 1, the eigenfunctions ^ i (u) are mapped into 

^ ( M ) = "lifrT^Xi ( u ) = V m ( u )Xi (u) , (A.IO) 

which explain the manifestation of the term \fm (u) in (41). Then due to the redef- 
inition of the scalar product, it is obvious that X \ ( u ) are the physical wavefunctions, 
while r i\) i (u) are called the group wavefunctions. 
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